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Weighted networks are more synchronizable: how and why 
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Most real-world networks display not only a heterogeneous distribution of degrees, but also a het- 
erogeneous distribution of weights in the strengths of the connections. Each of these heterogeneities 
alone has been shown to suppress synchronization in random networks of dynamical systems. Here 
we review our recent findings that complete synchronization is significantly enhanced and becomes 
independent of both distributions when the distribution of weights is suitably combined with the 
distribution of degrees. We also present new results addressing the optimality of our findings and 
extending our analysis to phase synchronization in networks of non-identical dynamical units. 
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I. INTRODUCTION 

The recent explosion in the study of complex networks 
of dynamical systems [l| has provided compeling evi- 
dence that the synchronization of coupled oscillators is 
drastically influenced by the structure of the underlying 

network of couplings 0,&SliH0ISH[l3- Prc " 
vious works on synchronization have focused mainly on 
the influence of the topology of the connections by assum- 
ing that the coupling strength is uniform. These works 
have shown that the ability to synchronize is generally 
enhanced in both small-world networks (SWNs) an( i 
random scale-free networks (SFNs) [l2| as compared to 
regular networks |l3j| . indicating that synchronizability 
is strongly related to the average distance between oscil- 
lators. More recently, it has been shown that synchro- 
nizability also depends critically on the heterogeneity of 
the degree distribution and that this can lead to counter- 
intuitive behaviors In particular, random networks 
with strong heterogeneity in the degree distribution, such 
as SFNs, are more difficult to synchronize than random 
homogeneous networks 0, despite the fact that hetero- 
geneity reduces the average distance between nodes . 

However, synchronization, as many other dynamical 
processes, is influenced not only by the topology, but 
also by the strength of the connections. Very recently, 
it has been shown that synchronization is suppressed 
in networks with heterogenous distribution of weights 
in the connection strengths Q, even when the distri- 
bution of degrees is homogeneous. This is an impor- 
tant result because most complex networks where syn- 
chronization is relevant arc indeed weighted and display 
a highly heterogeneous distribution of both degrees and 
weights 0,0,0- Examples include brain networks, 
both at the neuronal and cortical level [l^, an d airport 
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networks |l5j , which underlie the synchronization of epi- 
demic outbreaks in different cities Jljj . The identification 
and study of complex weighted networks with improved 
synchronization properties is thus of great interest. 

In this paper, we review this fundamental problem in 
the context of complete synchronization of identical oscil- 
lators |2fJ, |2l| . We introduce a weighted coupling scheme 
and we show that, for a given network topology, the syn- 
chronizability is maximum when the network of couplings 
is weighted and directed. For large sufficiently random 
networks, the maximum synchronizability is primarily 
determined by the mean degree of the network and does 
not depend on the degree distribution and system size. 
This contrasts with the case of unweighted (and undi- 
rected) coupling, where the synchronizability is strongly 
suppressed as the degree heterogeneity or number of os- 
cillators is increased. We also show that the total cost 
involved in the weighted couplings is significantly reduced 
as compared to the case of unweighted coupling and ap- 
pears to be minimum when the synchronizability is max- 
imum. Furthermore, we present new results that extend 
our findings to the synchronization of phase oscillators 
and provide evidence for the optimality of our weighted 
model. 

The paper is organized as follows. The problem of 
complete synchronization is formulated in Sec. HTland is 
analyzed in Sec. 11111 The problem of phase synchroniza- 
tion is considered in Sec. IfVI In Sec. we discuss the 
distribution of weights that would maximize synchroniz- 
ability for a given network topology. The conclusions are 
presented in the last section. 



II. COMPLETE SYNCHRONIZATION 

We introduce a weighted model of linearly coupled 
identical oscillators and we present a condition for the 
linear stability of the completely synchronized states in 
terms of the eigenvalues of the coupling matrix. 

The dynamics of a weighted network of N identical 
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oscillators is described by: 

N 

ii = F(x i )+a^^[H(x J -)-H(x i )], (1) 

N 

= F(x i )-<r^G ij H(x J -), i=l,...,^V, (2) 

i=i 

where F = F(x) governs the dynamics of each individual 
oscillator, H = H(x) is the output function, and a is 
the overall coupling strength. Matrix A — (Aij) is the 
adjacency matrix of the underlying network of couplings, 
where Aij = Wij if there is a link of strength wij > from 
node j to node i, and otherwise. Matrix G = (Gij), 

defined as Gij — —Aij + 5ij X) 7 =i Aiji 1S the coupling 
matrix. The rows of matrix G have zero sum and this 
ensures that the completely synchronized state {x;(i) = 
s(i), Vi | ds/dt = F(s)} is a solution of Eq. 

In the case of symmetrically coupled oscillators with 
uniform coupling strength, the network of couplings is 
unweighted and undirected, and G is the usual (symmet- 
ric) Laplacian matrix L — (Lij): the diagonal entries 
are La = fcj, where fcj is the degree of node i, and the 
off-diagonal entries are Lij = — 1 if nodes i and j are con- 
nected and = otherwise. For Gij — , heterogene- 
ity in the degree distribution suppresses synchronization 
in important classes of networks 7\. The synchronizabil- 
ity may be easily enhanced if we modify the topology of 
the network of couplings. Here, however, we address the 
problem of enhancement of synchronizability for a given 
network topology. 

In order to enhance the synchronizability of networks 
with heterogeneous degree distribution, we propose to 
scale the coupling strength by a function of the degree of 
the nodes. For specificity, we take 

Gij — Lij /k^ , (3) 

where (3 is a tunable parameter. We say that the network 
or coupling is weighted when (3 ^ and unweighted when 
(3 = 0. Networks with (3 = and (3 > are depicted in 
Figs.^a) and^Jb), respectively. 

The underlying network associated with the Laplacian 
matrix L is undirected and unweighted, but for (3 ^ 0, 
the network of couplings becomes not only weighted but 
also directed because the resulting matrix G is in general 
asymmetric. This is a special kind of directed network 
where the number of m-links is equal to the number of 
ow£-links in each node, and the directions are encoded in 
the strengths of in- and out-links [Fig. ^Jb)]. In spite of 
the possible asymmetry of matrix G, all the eigenvalues 
of matrix G are nonnegative reals and can be ordered as 
= Ai < A 2 • • • < Xn, as shown below. 

In matrix notation, Eq. J2Jl can be written as G = 
D^^L, where D = diag{fci, k%, . . . fejv} is the diagonal 
matrix of degrees. (We recall that the degree ki is the 
number of oscillators coupled to oscillator i.) From the 
identity det(D"' 3 L - XI) = det(L>-' 3 / 2 J LL>-' 3 / 2 - XI), 




FIG. 1: Illustration of model © for (a) j3 = (unweighted 
coupling) and (b) /3 > (weighted coupling). 



valid for any A, we have that the spectrum of eigenval- 
ues of matrix G is equal to the spectrum of a symmet- 
ric matrix defined as H — D~P/ 2 LD~P/ 2 . As a result, 
all the eigenvalues of matrix G are real, as anticipated 
above. Moreover, because H is positive semidefinite, all 
the eigenvalues are nonnegative and, because the rows of 
G have zero sum, the smallest eigenvalue Ai is always 
zero. If the network is connected, A2 > for any fi- 
nite (3. Naturally, the study of complete synchronization 
of the whole network only makes sense if the network 
is connected. For (3 = 1, matrix H is the normalized 
Laplacian matrix |22f . In this case, if N > 2 and the 
network is connected, then < A2 < N/(N — 1) and 
N/(N — 1) < Xn < 2. For spectral properties of un- 
weighted complex networks, see Refs. [22L l23l 124. l25| . 

The variational equations governing the linear stability 
of a synchronized state {x^(i) = s(i),Vi} of the system 
in Eqs. © and © can be diagonalized into N blocks of 
the form 

^ = [DF( S )-aDU(s)]r 1 , (4) 
where D denotes the Jacobian matrix, a = aXi, and 
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FIG. 2: Master stability function: illustration of the condition 
Qfl < 0X2 < .. . < CrAjV < Q2- 



Aj are the eigenvalues of the coupling matrix G. The 
largest Lyapunov exponent A(o) of this equation can be 
regarded as a master stability function, which determines 
the linear stability of the synchronized state [26|: the 
synchronized state is stable if A(aXi) < for i = 2, ... N. 
(The eigenvalue Ai corresponds to a mode parallel to the 
synchronization manifold.) 

For many oscillatory dynamical systems [|| |2(|, the 
master stability function A(o) is negative in a single, fi- 
nite interval (an, o 2 ). Therefore, the network is synchro- 
nizable for some a when ct\ < 0X2 < . . . < cf\n < 02, as 
depicted in Fig. [21 This is equivalent to the condition 

R = X N /X 2 < 0-2/0-1, (5) 

where 02/01 depends only on the dynamics (F, H, and 
s), while the eigenratio R depends only on the coupling 
matrix G. The problem of synchronization is then re- 
duced to the analysis of eigenvalues of the coupling ma- 
trix 0: the smaller the eigenratio R, the larger the syn- 
chronizability of the network and vice versa. Alterna- 
tively the eigenratio R can also be regarded as a measure 
of the speed to synchronize, where the speed is defined by 
the smallest transverse Lyapunov exponent. 



III. ENHANCEMENT OF 
SYNCHRONIZ ABILITY IN WEIGHTED 
NETWORKS 

We study the system in Eqs. J5J and © and we argue 
that the synchronizability is maximum (R is minimum) 
and the coupling cost is minimum for = 1. 



A. Mean Field Approximation 

A mean field approximation provides insight into the 
effects of degree heterogeneity and the dependence of R 
on 0. 

The dynamical equations (0 can be rewritten as: 

^=F(x i )+^ 1 - /3 [H i -H(xO], (6) 



where 

is the local mean field from all the nearest neighbors of 
oscillator i. If the network is sufficiently random and 
the system is close to the synchronized state s, we may 
assume that EL ~ H(s) and we may approximate Eq. © 
as 

^i=F(x l ) + ( rfc 4 1 -' 3 [H(s)-H(x 4 )], (8) 

indicating that all the oscillators are decoupled and 
forced by a common oscillator with output H(s). 

From a variational equation analogous to Eq. we 
have that all oscillators in Eq. (JHJ will be synchronized 
by the common forcing when 

Oi < ak\ 13 < o 2 Vi. (9) 

For ^ 1, it is enough to have a single node with degree 
very different from the others for this condition not to be 
satisfied for any a. In this case, the complete synchro- 
nization becomes impossible because the corresponding 
oscillator cannot be synchronized. This explains the re- 
sults of Ref. |3 on the suppression of synchronizability 
due to heterogeneity in unweighted networks. 

Within the mean field approximation, the eigenratio is 

R = (fcrnax/fcrnin) 1 "' 3 for < 1 and R = (U/fcmax) 1 ^ 

for > 1, where fc m in = nihiil&i} and A: max = maxijfci}, 
and reaches its minimum at = 1. This suggests that 
the maximum synchronizability in weighted networks is 
achieved when 0=1. As shown in Ref. 21], the same 
indication is provided by the analysis of a diffusion pro- 
cess related to the spread of information in the network. 
We thus expect the (un-approximated) system in Eq. @ 
to exhibit maximum synchronizability at = 1. 

We now test this prediction numerically on random 
SFNs |22j. The networks are generated as follows. Each 
node is assigned to have a number ki > fc m ; n of "half- 
links" according to the probability distribution P(k) ~ 
k^ 1 , where 7 is a scaling exponent and /c m i n is a constant 
integer. The network is generated by randomly connect- 
ing these half-links to form links, prohibiting self- and 
repeated links. We discard the networks that arc not 
connected. In the limit 7 — 00, all nodes have the same 
degree k = k mhl . 

Our numerical computations confirm that the eigenra- 
tio R has a well defined minimum at = 1, as shown in 
Fig.[3]for various 7. The only exception is the class of ho- 
mogeneous networks, where all the nodes have the same 
degree. When the network is homogeneous, the weights 
can be factored out and the eigenratio R is indepen- 
dent of (Fig. |3 solid line). Random homogeneous net- 
works correspond to random SFNs with 7 = 00. In SFNs, 
the heterogeneity increases as the scaling exponent 7 is 
reduced. As shown in Fig. |3 the minimum of the eigen- 
ratio R becomes more pronounced as the heterogeneity 
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FIG. 3: Eigenratio R as a function of j3 for random SFNs 
with 7 = 3 (•), 7 = 5 (■), 7 = 7 (a), and 7 = 00 (solid line). 
Each curve corresponds to an average over 50 realizations of 
the networks for fc m i n = 10 and N = 1024. 



FIG. 4: Eigenratio R as a function of the scaling exponent 7 
for random SFNs with (3 = 1 (•) and f3 = (o). The other 
curves are the approximation of the eigenratio in Eq. 1121 
(solid line) and the eigenratio at 7 = 00 (dashed line). The 
symbols correspond to random homogeneous networks with 
the same mean degree of the corresponding SFNs (the degrees 
are indicated in the figure). The other network parameters 
are the same as in Fig. [3] 



of the degree distribution is increased. A pronounced 
minimum for the eigenratio R at /3 = 1 is also observed 
in various other models of complex networks [20I l2lj . 



B. Mean Degree Approximation 

In what follows, we present an approximation for the 
eigenratio R. Here we focus on the case of maximum 
synchronizability {[3 — 1). 

Based on results of Ref. [24| for random networks with 
given expected degrees, we have 



c{l-A 2! A JV -l} = [1 



(10) 



Moreover, the semicircle law holds and the spectrum of 
matrix H is symmetric around 1 for fc min 3> Vk in the 
thermodynamic limit p4j . These results are rigorous for 
ensembles of networks with a given expected degree se- 
quence and sufficiently large minimum degree k mul , but 
our numerical computation supports the hypothesis that 
the approximate relations 
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hold under much milder conditions. In particular, rela- 
tions Ql 1[1 are expected to hold true for any large, suffi- 
ciently random network with fc m ; n 3> 1. 

Under the assumption that 1 — A2 ~ A^v — 1« 2/y/k, 
the eigenratio can be written as 



R 



1 + 2/Vk 
1 - 2/y/k' 



(12) 



Therefore, for /3 = 1, the eigenratio R is primarily deter- 
mined by the mean degree and does not depend on the 



number of oscillators and the details of the degree distri- 
bution. This is a remarkable result because, regardless of 
the degree distribution, the network at (3 = 1 is just as 
synchronizable as a random homogeneous network with 
the same mean degree, and random homogeneous net- 
works appear to be asymptotically optimal in the sense 
that R approaches the absolute lower bound in the ther- 
modynamic limit for large enough k 0. 

We now test our predictions on random SFNs. As 
shown in Fig. in unweighted SFNs, the eigenratio R 
increases with increasingheterogeneity of the degree dis- 
tribution (see also Ref. Q). But, as shown in the same 
figure, the eigenratio does not increase with heterogene- 
ity when the coupling is weighted at /3 = 1. The differ- 
ence is particularly large for small scaling exponent 7, 
where the variance of the degree distribution is large and 
the network is highly heterogeneous (note the logarith- 
mic scale in Fig. 0J. For (3 = 1, the eigenratio R is well 
approximated by the relation in Eq. (|12J) [Fig. 0] solid 
line]. For (3=1, the eigenratio of the SFNs is also very 
well approximated by the eigenratio of random homoge- 
neous networks with the same number of links [Fig. 
0]. Therefore, for (3 = 1, the variation of the eigenratio 
R with the heterogeneity of the degree distribution in 
SFNs is mainly due to the variation of the mean degree 
of the networks, which increases in random SFNs as the 
scaling exponent 7 is reduced. 

In Fig. El we show the eigenratio R as a function of the 
system size N . In unweighted SFNs, the eigenratio in- 
creases strongly as the number of oscillators is increased. 
Therefore, it may be very difficult or even impossible 
to synchronize large unweighted networks. However, for 
[3 = 1 , the eigenratio of large networks appears to be in- 
dependent of the system size, as shown in Fig. [5] for ran- 
dom SFNs. Similar results are observed in many other 
models of complex networks. Altogether, these provide 
strong evidence for relation l|12|) and shows that synchro- 
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FIG. 5: Eigenratio R as a function of the number of oscillators 
for random SFNs with 7 = 3. Dotted lines are guides for the 
eyes. The legend and other parameters are the same as in 
Fig-H 
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FIG. 6: Normalized cost C/Nai as a function of f3 for random 
SFNs with scaling exponent 7 = 3. The other parameters are 
the same as in Fig. El 



FIG. 7: Normalized cost C/Nai as a function of the scaling 
exponent 7 for random SFNs with /3 = 1 (•) and j3 = (o), 
and for random homogeneous networks with the same mean 
degree (0). The dashed line corresponds to 7 = 00. The other 
parameters are the same as in Fig. 13 



approximated by the cost for random homogeneous net- 
works with the same mean degree [Fig.0 0] , in agreement 
with our analysis in Sec. lIIIBl that. at = 1, the eigen- 
value A2 is fairly independent of the degree distribution. 
The cost at j3 = 1 is significantly reduced as compared 
to the case of unweighted coupling (f3 = 0), as shown in 
Fig. The difference becomes more pronounced when 
the scaling exponent 7 is reduced and the degree distri- 
bution becomes more heterogeneous. Similar results are 
observed in other models of complex networks. There- 
fore, cost reduction is another important advantage of 
suitably weighted networks. 



nizability is significantly enhanced for /? = 1 as compared 
to the case of unweighted coupling ((3 = 0). 



C. Coupling Cost 

We now turn to the problem of the cost involved in 
the network of couplings. The total cost C is naturally 
defined as the total input strength of the connections of 
all nodes at the synchronization threshold: 

N 

C = a min J2kl~^, (13) 

i=i 

where er m i n = ai/A2 is the minimum coupling strength 
for the network to synchronize. We recall that a\ is the 
point where the master stability function first becomes 
negative. For j3 = 1, we have C = Na\/\2- 

As a function of (3, the cost has a broad minimum at 
(3 = 1, as shown in Fig. for random SFNs. This result 
is important because it shows that maximum synchro- 
nizability and minimum cost occur exactly at the same 
point. The cost for random SFNs at j3 = 1 is very well 



IV. PHASE SYNCHRONIZATION 

The above analysis on identical oscillators can serve 
as a good approximation even when the oscillators are 
not fully identical. In some realistic situations j28J , how- 
ever, the interacting oscillators are significantly different, 
specially in the oscillation frequencies |2^. Here we ex- 
tend the enhancement of synchronizability to networks 
of non-identical oscillators. 

We consider phase synchronization in complex net- 
works of limit cycle oscillators, 

N 

i i = u i F(xi)-a^G ij U{x j ), i = l,...,N, (14) 
i=i 

where x = F(x) represents a Van der Pol oscillator, u = 
v, v = 0.3(1 — u 2 )v — u, and the couplings are through 
both variables u and v such that H(x) = x. We assume 
that the frequencies LOi follow a uniform distribution in 
an interval [1 — Acj, I + Aw]. In our simulations we set 
Aui = 0.2. 

Now we study the effects of the weighted coupling 
on phase synchronization. For an arbitrary /3, on aver- 
age each oscillator receives an effective input of strength 
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FIG. 8: Time series of the mean field X in random SFNs 
of non-identical Van der Pol oscillators with 7 = 3. The 
coupling strengths are (a) a* — 0.1 and (b) a* = 0.2. The 
other parameters are the same as in Fig. 
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FIG. 9: Amplitude as a function of a* for random SFNs of 
Van der Pol oscillators with (3 = (o) and (3 = 1 (•). The 
symbols correspond to random homogeneous networks with 
the same mean degree of the random SFNs. The results are 
averaged over 20 realizations. The scaling exponent is 7 = 3 
and the other parameters are the same as in Fig. |31 



V. WHAT ARE THE MOST 
SYNCHRONIZABLE NETWORKS? 



a* = a k\ ~^/N. We can compare the synchroniza- 
tion behavior at different (3 for the same effective strength 
a*. In Fig. |H1 we show time series of the mean field 
X = Yli=i x i/N for unweighted ((3 = 0) and weighted 
random SFNs ((3 = 1). For small coupling strength a*, 
neither of the networks display significant collective be- 
havior [Fig.|Sfa)]. As the coupling strength is increased, 
collective oscillations emerge for both weighted and un- 
weighted networks, but the oscillations are much more 
pronounced for the networks with (3=1 [Fig. |Sfb)]. 

In Fig. we show the amplitude of the mean field 
X as a function of a*. Here the amplitude is defined 
as the standard deviation of X over time. The ampli- 
tude is approximately zero for small coupling strength, 
increases sharply as a* is increased beyond a certain crit- 
ical value, and saturates for large a* (Fig.|5J). (This tran- 
sition becomes sharper as the number N of oscillators is 
increased). The overall behavior is similar for weighted 
and unweighted networks but, again, the amplitude is 
significantly larger for networks with = 1 at the same 
value of a*. Moreover, the amplitude of random SFNs 
at j3 = 1 is well approximated by the amplitude of ran- 
dom homogeneous networks with the same mean degree 
(Fig. [HJ) , which are networks that exhibit good phase syn- 
chronization properties. 

All these provide strong evidence that our weighted 
coupling scheme strongly enhances the synchronizability 
of networks of non-identical oscillators. This is expected 
to be relevant for realistic networks. 



Within the weighted coupling scheme of Eq. ©, the 
networks that exhibit maximum synchronizability and 
minimum cost are those with (3 = 1 fSec. IIII(I . Physically, 
this seems to be related to the fact that the synchroniza- 
tion of the network is determined by the input signal at 
each node (as opposed to the output signal) and that the 
sum J^j Aij = fcl ^ of the strengths of all the in- links 
of a node becomes uniformly equal to 1 for all the nodes 
exactly at (3 = 1. This leads to the natural question of 
whether there is any other distribution of weights that 
can further improve the synchronizability and cost of the 
network for a given topology of connections. As we show, 
the answer is positive for finite size networks, but seems 
to be negative in the thermodynamic limit. 

In order to address this question we introduce the fol- 
lowing weighted coupling scheme. If there is a link be- 
tween a node i with degree ki and a node j with degree 
kj , the strength Aij = Wij of link j — * i is taken to be 

Wtj = (hkjf/Nr (15) 

where 8 is a tunable parameter and N r = ^Zj^ihikj) is 
a normalization factor that ensures that Y] j Aij = 1 (the 
sum i s over an the neighbors of node i). For 8 = 

we recover the model in Eq. © for (3=1. 

This new weighted model is partially motivated by 
observations in real networks, including scientific col- 
laboration networks Il5l , airport networks |15l Il6| , and 
metabolic networks 16], which have been found to ex- 
hibit a strong correlation between the weight Wij and 
product of the corresponding degrees as (wy) ~ (kikj) 8 . 
A similar model has been considered in the study of co- 
herence in networks of chaotic maps Q. 
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FIG. 10: (a) Eigenratio R and (b) normalized cost C/{Not\) 
as lunctions ol 9 for random SFNs with with 7 = 3 (•), 7 = 5 
(■), and 7 = 7 (a). The other parameters are the same as in 
Fig. 




FIG. 11: The same as in Fig^for random SFNs with 7 = 3 
and various choices for the number N of oscillators, as indi- 
cated in the figure. 



We now study the weighted model l|15|) numerically 
in the context of complete synchronization of identical 
oscillators. In Fig. ^3 we show the eigenratio R and 
the normalized cost C/{Not\) as functions of 9 for dif- 
ferent values of the scaling exponent. In each case, both 
the eigenratio and the cost reach a minimum for some 
6*(N, r y) > 0, which is approximately the same for R 
and C and indicates that these networks synchronize 
better than those of the weighted model |J3} for /3 = 1 
(9 = 0). However, as compared to the difference between 
unweighted networks (/3 = 0) and networks weighted for 



(3 = 1 in model (J3J), the further reduction of R and C 
within model l|15l) is very small. 

Indeed, R(J3 = O)-R(0 = 1) and C(/3 = 0)-C(/3 = 1) 
[Fig. 0] are orders of magnitude larger than 
R(0 = 0) - R(8 = 9*) and C{9 = 0) - C(9 = 8*) 
[Fig. even when the networks are strongly heteroge- 
neous (i.e., 7 is small). Moreover, the difference between 
the synchronizability at 9 — and 9 — 9* reduces when 
the number N of oscillators is increased, as shown in 
Fig. ^2 Since 9* approaches 9 = for increasing N, 
our numerics are consistent with the hypothesis that the 
minimum of R and C will be at 9 = for N — > 00. That 
is, the weighted model (0) for (5 = 1 may be optimal 
in the thermodynamic limit. The verification of this 
hypothesis is currently an interesting open problem. 

VI. CONCLUSIONS 

Complex networks with strong heterogeneity in the dis- 
tributions of either degrees or weights are poorly synchro- 
nizable. This suppression of synchronizability is due to 
the different intensities of the input signals received by 
different nodes in the network. We have shown that the 
synchronizability is significantly improved when the dis- 
tribution of weights is constrained by the distribution of 
degrees in such a way that all the nodes receive the same 
intensity of input signal. In this case, the distributions of 
degrees and weights as well as the distribution of the to- 
tal strength of out-links per node, ^ i A+j, may be highly 
heterogenous, but the distribution of the total strength 
of in-links, y\ Aij, is the same for all the nodes in the 
network. In particular, our analysis shows that properly 
weighted complex networks are much more synchroniz- 
ablc than the corresponding unweighted networks. 
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